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ABSTRACT 

We calculate the minimum value of the power in kinetic bulk motion of the galactic 
superluminal source GRS 1915+105. This value far exceeds the Eddington luminosity 
for accretion onto a black hole of 10 solar masses. This large value severely limits the 
possible carriers of the kinetic luminosity at the base of the jet, and favours a jet 
production and acceleration controlled by a magnetic field whose value, at the base of 
the jet, exceeds 10 8 Gauss. The Blandford and Znajek process can be responsible of 
the extraction of the rotational energy of a Kerr black hole, if lasting long enough to 
provide the required kinetic energy. This time, of the order of a day, implies that the 
process must operate in a stationary, not impulsive, mode. 
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1 INTRODUCTION 

GRS 1915+105 was discovered in 1992 with the WATCH 
telescope on board the GRANAT satellite (Castro-Tirado 
et al. 1992). A radio counterpart was subsequently iden- 
tified and bipolar outflows with apparent superluminal mo- 
tions were observed; the standard interpretation of this phe- 
nomenon in terms of relativistic jets (Rees 1966) places the 
source at a distance D=12.5 kpc at an angle i — 70° to 
the line of sight (Mirabel & Rodriguez 1994). The source 
cannot be observed in the optical band, due to the heavy 
extinction in the Galactic plane, but the X-ray luminosity 
often well above the Eddington luminosity for a neutron 
star suggests the presence of a ~ 10 solar mass black hole. 
Other circumstantial evidence comes with the similarities 
with the other galactic superluminal source GRO J1655-40, 
which has been shown unambiguously to harbor a compact 
object of ~ 7Mq (Orosz & Bailyn 1997). Since its discov- 
ery, GRS 1915+105 has displayed an extraordinary richness 
in variability in the X-ray band (Greiner, Morgan & Remil- 
lard 1996, Chen, Swank & Taam 1997, Belloni et al. 1997), 
in the IR and radio bands (Fender et al. 1997, Pooley & 
Fender 1997). Recent (Eikenberry et al. 1997, Mirabel et al. 
1997) simultaneous multiwavelength observations of this su- 
perluminal source show a strict link between activity in the 
inner accretion disk and plasma ejections. 

In this paper we estimate the minimum kinetic power 
associated with the ejections and investigate the possible 
energy transport mechanism from the compact source to 
distances where the radio blobs are resolved, showing that 
the more suitable solution implies the presence of a strong 
magnetic field at the footpoint of the jet. 



In section 2 we summarize the existing radio and IR 
data of the ejection events. In section 3 we directly calculate 
the minimum kinetic power during the ejection. In section 
4 we investigate the possible origin of the kinetic power. In 
section 5 we summarize the main results of the paper and 
discuss their consequences. 



2 OBSERVATIONAL CONSTRAINTS DURING 
THE OUTFLOWS 

Radio observations show that GRS 1915+105 in its active 
state is characterized by the increase of the flux level from 
S v < 10 mjy to S v ~ 100 mjy, the so called plateau state, 
on the top of which radio flares are superimposed, with du- 
rations from days to a month and fluxes up to a few Jy 
(Foster et al. 1996). The spectra during the plateau state 
are generally flat or inverted, a < (with S v oc v~~ a ), prob- 
ably indicating synchrotron self-absorption. On the other 
hand the spectra of radio flares reveal a transition from op- 
tically thick to thin (a > 0) emission during the rise stage, 
with spectra harder (a ~ 0.5) when the fluxes reach the 
maximum and softer (a ~ 1) during the subsequent decay, 
generally interpreted as indicator of synchrotron radiation 
of an expanding radio cloud. 

Observations with the Very Large Array and multiwave- 
length campaigns (with VLA and UKIRT simultaneously) 
pointed out that some strong radio flares are related with 
the ejection of radio clouds. There are two different kind 
of ejections: 1) major ejections as the prominent radio out- 
burst of March 1994 (Mirabel & Rodriguez 1994), where 
the ejecta actually moved for several weeks along a direc- 
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tion forming an angle 6 — 70° to the line of sight, with 
bulk velocity of 0.92c and expansion at ~ 0.2c, with a spec- 
tral index change from a — 0.49 (when the blobs could not 
be resolved) to a = 0.84 (when the two condensations had 
moved apart); and 2) so-called "baby jets" (Eikenberry et al. 
1997, Mirabel et al. 1997): quasi-periodic oscillations in the 
radio flux, coupled with similar oscillations in X-ray and IR 
bands and interpreted as scaled-down (in space and time) 
ejections of synchrotron emitting expanding blobs. The anal- 
ysis of QPOs in IR and radio bands clearly shows that there 
is a time shift of flux peaks, with short wavelengths peak- 
ing first: the IR (2 am) peak precedes the 2 cm peak, that 
precedes the 3.6 and 6 cm peaks. 

It must be emphasized that in one occasion (July 1995) 
also in the near-infrared K band (A = 2.2/im) a jet was 
observed, with the same position angle of the radio jet, and 
near-infrared magnitude of A"=13.9 (slightly brighter than 
the total source magnitude at its weakest measured value 
of K=1A.Z) and separated from the central source by 0.3" 
(Sams et al. 1996). 



3 MINIMUM KINETIC POWER CONDITION 

In order to deduce the minimum kinetic luminosity Lk re- 
lated to a major ejection we consider the March 1994 event 
for which the observational data are the most exhaustive. 
Since we are interested more in power rather than in energy 
estimates, we follow an alternative way with respect to the 
standard one, that consists to determine first the internal 
energy of the blob via the minimum energy criterium and 
then to obtain an estimate of the kinetic power, by making 
some (highly uncertain) assumptions about the energization 
time. 

We instead directly calculate the kinetic power, equal 
to the energy flux through a cross section of the jet (see 
also Ghisellini & Celotti, 1998). This energy flux can be 
carried by particles and by the toroidal magnetic field, and 
correspondingly we have 



Lk, P = nR 2 T 2 Pen < 7 > m e c 2 ( 1 H ^± I 

\ m e < 7 > J 



L k , B = hR 2 T 2 I3cUb 
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where m+ is either m p in the case of a "normal" e— p plasma 
or < 7 > m e + for e^ pairs, T is the Lorentz factor of the 
bulk motion, R is the cross section radius of the jet, n' in 
the comoving particle density, Ub is the magnetic energy 
density measured in the comoving frame and < 7 > is the 
mean Lorentz factor of the electrons, as measured in the 
comoving frame. A lower limit on n' can be estimated by the 
observed synchrotron emission L syn of the blob. Assuming 
a spherical emitting volume of radius R, the number density 
of leptons producing the observed radiation is 
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where < j 2 > is averaged over the relativistic electron dis- 
tribution, and 8 = [r(l — /3cos^)] _1 is the Doppler factor. 
The B~ 2 dependence of the estimated particle density al- 
lows to minimize the total power Lk, P + Lk,B with respect 



to the magnetic field, since both the bulk Lorentz factor F 
(—2.55) and the viewing angle 9 (=72°) are known: 
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yielding the value of the magnetic field B„ 
to the minimum power: 
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where / is a parameter directly related to the shape of the 
electron energy distribution, and therefore to the shape of 
the synchrotron radiation they emit: 
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Since the radio emission is a power law with spec- 
tral index a ~ 0.5, the particle distribution N(-y) oc 
7~ 2 between 71 and 72, yielding / ~ ^(72/71 )/72{l + 
m+/m e [7i ^(72/71)]}. With R ~ 7 x 10 15 cm (geometric 
average of the size of blob of GRS 1915+105, assuming a 
distance of 12.5 kpc), L eyn = 10 33 erg s _1 and j3 = 0.92 



(corresponding to V = 2.55), we obtain B v 



0.12 G 



and B min>e ± = 0.036 G if / = 1.8 (71 = 1 and 72 = 10 3 , as 
required for production of the observed synchrotron photons 
at v ~ 300 GHz). A low energy cut-off in the particle dis- 



tribution would decrease somewhat B„ 



= 0.054 



G and B min e ± = 0.031 G for 71 = 30), while the high en- 
ergy cut-off is consistent with the production of the observed 
radiation at 300 GHz (with B = B min ). However, the depen- 
dence of B m in on the extremes of the electron distribution 
is rather weak (B m in oc f 1 ^ 4 ). 

Note that the value B m in,e P = 0.12 is almost a factor 
3 greater than the one obtained by Mirabel & Rodriguez 
(1995) and Liang & Li (1995), while it substantially agrees 
with the value estimated from the requirement of optical 
transparency of the plasmoid with respect to synchrotron 
self absorption (Atoyan & Aharonian 1997). 

With a magnetic field value equal to B m in , the particle 
kinetic power and the Poynting flux are nearly equal, and 
the total power 



3r 2 /3c 
28 2 



irni e c 2 RL syn f 



arc 



1/2 



(6) 



is, respectively for e — p plasma and e ± pairs, 

~ 3.3 x lO^ergs" 1 and ~ 2.9 x 10 39 ergs _1 

This is the minimum kinetic luminosity involved in major 
ejection events, calculated assuming 71 = 1 and 72 = 10 3 



(Lk,tot,e 



1.4 x 10 3a erg s _i and L k tot e ± ~ 2.1 x 10 3a erg 



s- 1 if 7l =30). 

Setting B — Bmin in equation (2) the total numbers 
of emitting particles in the magnetic cloud are respectively 
N' e _ = 5.9 x 10 47 and N' e± = 6.9 x 10 48 , corresponding to 
internal energies 



7- 



Ei = rn e c 2 / 7_/V(7)d7 = N < 7 > m e c 2 



~ 3.3 x 10 42 erg, 



3.9 x 10 43 erg 



(7) 
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These internal energies can be compared with the corre- 
sponding kinetic energies in bulk motion E k given by 

R 



Ek = Lkt ~ Lk — 



8.5 x 10 45 erg 



and 



7.4 x 10 44 erg (8) 



where t is the time needed for the blob to cross the jet 
section. 



Compton cooling time and the dynamical time Rd/c for a 
particle of Lorentz factor 7 is 

R d /c cr T L d 7 

A large fraction of Lk would be lost and converted into ra- 
diation, contrary to what observed. We therefore conclude 
that a pure e^ 1 pair plasma cannot carry the kinetic power, 
irrespective of whether it is hot or cold. 



4 ORIGIN OF THE KINETIC LUMINOSITY 

The kinetic power calculated above refers to the radio emit- 
ting blob, at a huge distance from the putative black hole 
and accretion disk. The most economic assumption is that 
this power is approximately conserved along the jet, with 
very small dissipation giving rise, e.g., to the random par- 
ticle energy responsible for the emission. If not, we would 
have to assume that at the base of the jet an even larger 
value of Lk is produced. In principle there are several possi- 
ble energy carriers: e ± pairs, p — e~ , pure magnetic field or 
a mixture of these components. Consider also that the par- 
ticles can in principle be "hot" or "cold" : at large distances 
at least some of the electrons are "hot" (indeed, in our es- 
timates we assumed that all the electrons are "hot", since 
this is the most economic assumption), however these "hot" 
particles cannot come directly from the inner region (see be- 
low), but need to be accelerated or reaccelerated, therefore 
for the inner jet energy carriers one has to consider both 
cases. Furthermore, we have to note that in the "hot" case, 
the particle random energy, increasing the mass, can con- 
tribute to the kinetic power. Let us examine the possible 
cases. 



4.1 "Cold" e ± pairs 

If the kinetic power is carried by a pure e ± cold pair plasma, 
we can calculate the corresponding pair density and scatter- 
ing optical depth at some jet radius R close to the base of 
the jet 

crLk 



t± = arRn± — 



7ri?r 2 /3m e c 3 



4.2 x 10 



Lk 



10 7 cm 



2.9 x 10 39 erg s" 1 R ^ 

since the annihilation time-scale is of the order of R/(cr±), 
cold pairs cannot survive annihilation. 



4.2 "Hot" e ± pairs 

If the pairs are relativistic, with average random Lorentz 
factor < 7 >, the above estimate of t± decreases by a factor 
< 7 >. In addition, the annihilation cross section decreases. 
However these relativistic pairs are embedded in a dense 
radiation field, produced by the accretion disk, and they 
cool on time-scales shorter than the dynamical time-scale. 
In a region of size R d the radiation energy density due to the 
accretion disk luminosity Ld produced within Rd is of the 
order of Ud = L d /{^R 2 d c). The ratio between the inverse 



(11) 



4.3 "Normal" e — p plasma 

The case of electrons with average random Lorenz factor 
< 7 > greater than m p /m e and cold protons is analogous to 
the above case: inverse Compton scattering of seed accretion 
disk photons cools the electrons in a time shorter than Rd/c. 

The case of "hot", relativistic protons is instead im- 
mune to radiative losses, and the main concern comes from 
the needed confining pressure. If the kinetic power is car- 
ried by protons with an average Lorentz factor < -y p >, the 
corresponding (comoving) pressure is of the order of their 
energy density. If the confinement is magnetic, the required 
value of the magnetic field is 

B ( 8Lk 

\R 2 r 2 (] C/ 

which is equal to the value that the magnetic field must have 
if it is the main carrier of the kinetic power. 

If both electrons and protons are cold, the electron scat- 
tering optical depth is a factor m p /m e smaller than the one 
of equation (9), i.e. of the order of a few. There are no severe 
limits in this case: the Compton drag is not sufficient to slow 
down the plasma in the jet. It is however useful to estimate 
the predicted power emitted by the Compton drag process, 
because it can turn out to be important in sources where 
the jet is more aligned with the line of sight. For this simple 
estimate we will assume that the jet is Compton thick (i.e. 
tt > 1), so that the effective cross section of the process is 
the geometrical one, i.e. ttR 2 . As before, we assume that the 
accretion disk produces the luminosity L d within the region 
i?d > R, at the typical frequency v d - With these hypotheses, 
the observed bulk Compton luminosity Lbc is 



Lie 



(-Rd) 



Ld^2 c4 



(12) 



Most of this luminosity is observed at the frequency v ~ 
STv d . With R/R d ~ 1/10, L d ~ 10 39 erg s" 1 and <5 = 0.57 
we obtain Ltc ~ 1-7 x 10 erg s _1 and v ~ v d . This ra- 
diation is therefore unobservable in the case of the known 
galactic superluminals, which are both observed with large 
viewing angles, but may be very important for sources with 
5 ~ T > 1 (see also Sikora et al., 1997). 

4.4 Poynting vector 

The value of the magnetic field needed to carry Lk in the 
vicinity of the black hole is equal to equation (11), if the 
bulk of the magnetic field is moving with T, while the B- 
value before acceleration is a factor T greater, corresponding 
to B ~ 3 x 10 8 G at R = 10 7 cm. 

It must be remarked that the standard theory of ac- 
cretion disks (Shakura & Sunyaev 1973) predicts that the 
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maximum possible magnetic field at a given R/Rs (where 
Rs is the Schwarzshild radius) in a radiation dominated disk 
is B oc (M/Mq)~ 1//2 , so that microquasars will have mag- 
netic fields ~ 10 4 stronger than those in quasars. Therefore, 
the theoretical estimate B ~ 10 2 — 10 4 G obtained with the 
Blandford & Payne model (Blandford & Payne 1982) near 
the massive black holes in active galactic nuclei, clearly show 
that in microquasars B ~ 10 8 G can easily be attained. 

Note that this value of the magnetic field would be of 
the same magnitude of the magnetic field required by the 
Blandford-Znajek (1977) process to produce L k by extract- 
ing the rotational energy of a 10 solar mass Kerr black hole: 

L rot = 10 41 B%M\ erg s" 1 (13) 

where B — 10 9 i?9 G and M — 10Mi solar masses. 

From the above arguments we conclude that e or rel- 
ativists electrons with < 7 > greater than m p /m e cannot 
carry Lk in the inner region of the jet, while protons and 
the Poynting vector can. 

An alternative solution could be that the energy, carried 
in the inner region by "normal" cold plasma or toroidal mag- 
netic field, is converted in e ± pairs at large distance from the 
black hole and accretion disk. Pairs are produced most effi- 
ciently through photon-photon collisions, but this requires 
a powerful 7-ray continuum. Since the maximum efficiency 
in converting pairs in this way is of the order of 10 per cent 
(Svensson 1987), the high energy luminosity has to exceed 
10 40 erg s _1 , to produce a kinetic power of the order of 10 39 
erg s _1 . This luminosity is not observed. Therefore we con- 
clude that e ± pairs do not play any role as energy carriers 
along the jet and the minimun kinetic luminosity involved 
in major ejection events is given by the value estimated for 
e — p plasma (~ 3 x 10 40 erg s^ 1 ). 

There are no strong observational arguments to decide 
between protons and the Poynting vector, unless we observe 
a new galactic superluminal at a small viewing angle. Note 
however that if Lk is initially carried by a large magnetic 
field there is the possibility to tap a great reservoir of energy 
(the rotational energy of the black hole) and to accelerate 
the plasma to relativistic speeds. 

We therefore conclude that the kinetic power carried 
initially by the magnetic field is the most economic way to 
explain the observed energetics. Differential accretion disk 
rotation or rotating black holes can amplify magnetic fields 
to the required values, and the magnetic field will then ac- 
celerate particles to relativistic speeds. 

The amount of matter present at the base of the jet 
may not be negligible, especially if magnetic field lines help 
to channel particles from the disk to the jet. We can compare 
the inflow rate Mi n of the accretion process with the out- 
flow rate M out necessary to account for the kinetic power, in 
particles, of the radio blob. If we write the accretion lumi- 
nosity, as usually, as L acc — riM in c 2 , and the kinetic power 
as Lk = (r — l)M out c 2 , we derive 

Mont = M in (14) 

1 1 J^acc 

Since r\ ~ 0.1 and Lk ~ 10L acc we have that the the inflow- 
ing and the outflowing mass rates are comparable. This in 



turn suggests that most of the matter in the jet may come 
from the accretion disk. 



5 THE EJECTION TIME 

Another important question concerns the duration of the 
ejection events. When the blob becomes visible in the radio, 
it has a size of a few light days. It seems unlikely that it 
corresponds to an ejection duration lasting more than that, 
while a shorter ejection time may be possible. In the latter 
case the requirement on the initial kinetic power correspond- 
ingly increases. Note that direct radio observations of flare 
events do not directly constrain the ejection time, since the 
radio flux eventually produced in the first parts of the jet is 
heavily self absorbed. Again, the minimum power require- 
ments criterium favours ejection events which last for t out ~ 
few days. This corresponds to ~ 10 9 Rs,m/c, where Rs,io 
is the Schwarzchild radius for a 10 solar mass black hole. 
The immediate consequence is that the ejection event must 
be considered as a continuous and stationary process. Scaling 
for a superluminal active galactic nucleus, we would have, 
for a 10 9 solar mass black hole, an ejection phase lasting for 
2 x 10 5 years. 



6 DISCUSSION 

We have obtained a reliable lower limit to the kinetic 
power corresponding to major ejection events in superlumi- 
nal galactic sources, in particular for GRS 1915+105. This 
limit is of the order of 3 x 10 40 erg s _1 , much greater than the 
observed radiative luminosity, which is probably Eddington 
limited to values of the order of 10 39 erg s _1 , corresponding 
to a black hole of 10 solar masses. This by itself suggests 
that the jet acceleration mechanism cannot be radiative. 

We have investigated the role of e pairs, of normal 
plasma and of the magnetic field as energy carriers of the 
kinetic power in the inner jet regions, excluding an impor- 
tant role for the e ± pairs, and favouring a scenario in which 
the inital acceleration phase is controlled by a magnetic field 
of the order of 10 Gauss or more, able to channel and ac- 
celerate accretion disk matter in the jet. 

Conservation of kinetic power dictates that the injection 
time-scale is of the order of a day, a very long time-scale if 
measured in units of the light crossing time of a 10 solar mass 
Schwarzchild radius, indicating a stationary, not impulsive, 
process. Shorter injection times, albeit possible, correspond 
to larger kinetic powers, exacerbating the problem of how to 
obtain them. The rough equality between the value of the 
magnetic field needed to carry the kinetic power in the in- 
ner jet and the value necessary to tap the rotational energy 
of a Kerr black hole via the Blandford & Znajek process 
can be regarded as circumstantial evidence that this process 
is indeed the one responsible for the jet formation and ac- 
celeration. This process is candidate to power also the jets 
in radio-loud quasars, and we can compare the results ob- 
tained here with the corresponding estimates of the kinetic 
power of AGN, as derived by Celotti, Padovani & Ghisellini 
(1997) for a sample of radio-loud sources. For those AGNs, 
the above authors find a kinetic power between 10 45 and 
10 48 erg s _1 , of the same order of the luminosity needed 
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to ionize the broad line region of the same objects, and in 
agreement with the power required by the existence of the 
outer radio lobes. In the AGN case the kinetic and the ac- 
cretion luminosity are roughly equal, while in the galactic 
superluminal objects the kinetic power dominates. Another, 
obvious, difference is the Lorentz factor of the bulk motion, 
much greater in the AGN case. This results in a ratio of the 
outflowing to infalling mass rate M ou t/Mi n of order unity 
for the galactic superluminals and two orders of magnitude 
smaller for AGN. These estimates makes the galactic super- 
luminal sources the most efficient engines to produce colli- 
mated relativistic bulk motion, with the possible exception 
of gamma-ray bursts. 
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